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Abstract 

The nonlinear Schrodinger equation has several families of quasi-periodic travelling waves, 
each of which can be parametrized up to symmetries by two real numbers: the period of the 
modulus of the wave profile, and the variation of its phase over a period (Floquet exponent). In 
the defocusing case, we show that these travelling waves are orbitally stable within the class of 
solutions having the same period and the same Floquet exponent. This generalizes a previous 
work 22 where only small amplitude solutions were considered. A similar result is obtained 
in the focusing case, under a non-degeneracy condition which can be checked numerically. The 
proof relies on the general approach to orbital stability as developed by Grillakis, Shatah, and 
Strauss ^1 E], and requires a detailed analysis of the Hamiltonian system satisfied by the 
wave profile. 



Running head: Periodic waves in the NLS equation 

Corresponding author: Thierry Gallay, Thierry.Gallay@ujf-grenoble.fr 

Keywords: Nonlinear Schrodinger equation, periodic waves, orbital stability 



1 Introduction 



This paper is devoted to the stabihty analysis of the quasi-periodic travelhng wave solutions of the 
cubic nonlinear Schrodinger (NLS) equation 

iUt{x,t) + Ua:x{x,t) +j\U{x,t)fU{x,t) = , x^R, teR, (1.1) 

where 7 G {—1; 1} and U{x,t) £ C. Ea. Hl.lf) is a universal envelope equation describing the prop- 
agation of weakly nonlinear waves in dispersive media (see [20] for a comprehensive introduction). 
The nonlinearity in (jLlj) is "attractive" if 7 = -|-1 (focusing case) and "repulsive" if 7 = — 1 
(defocusing case). In both cases Eq. (|l.H) has a family of quasi-periodic travelling waves of the 
form 

U{x, t) = e'^P^-"^*) V{x-ct) , xeR, t € M , (1.2) 

where p,uj,c are real parameters and y : M ^ C is a periodic function. The simplest elements of 
this family are the plane waves, for which V is identically constant and p, oj satisfy the dispersion 
relation = w -|- 7|V^p. It is well-known and easy to verify that the plane waves are dynamically 
stable in the defocusing case, and unstable (if ^ / 0) in the focusing case [HI El- We shall 
therefore concentrate on the less explored situation where y is a nontrivial periodic function. In 
that case, we shall refer to ()1.2|) as a periodic wave, although U{x, t) is in general a quasi-periodic 
function of both x and t. 

The number of parameters in H1.2() can be reduced if we use the symmetries of Ea. Hl.ll) . We 
recall that the NLS equation is invariant under the following transformations: 

(i) U{x,t) I— > U{x,t)e^'^, if £ R (Phase invariance); 

(ii) U{x,t) ^ U{x -|- ^,i), ^ G M (Translation invariance); 

(iii) U{x,t) ^ e-'^^''+^^^U{x + vt,t), V gR (Galilean invariance); 

(iv) U{x,t) 1-^ \U{\x, X^t), A > (Dilation invariance). 

If U{x,t) is a periodic wave as in H1.2|) . we can use the Galilean invariance to transform it into a 
solution of the same form with c = 0. Then, using the dilation invariance, we can further assume 
that uj £ {-1;0; 1}. It follows that U{x,t) = e"''^*Ty(j;), where W{x) = e'PW{x) is a solution of 
the ordinary differential equation 

W^^{x) + ioW{x) + -f\W{x)\^W{x) = , xeR. (1.3) 

The bounded solutions of (\l.'d\i are completely classified for all values of the parameters uj,^. 
The simplest ones are the plane waves W{x) = Ae^^^ where p E M, A E C and p'^ = lo + "y\A\'^. 
The periodic waves correspond to quasi-periodic solutions of (|1.3|) of the form PF(x) = r{x) e^'^^^\ 
where r, ip are real functions with the property that r and <px are periodic with the same period. 
It turns out that Eq. (|1.3|l has a four-parameter family of such solutions, both in the focusing and 
in the defocusing case (see Sections |21 and [l] below). Actually, if 7 = —1, we must assume that 
id = 1 otherwise (|1.3j) has no nontrivial bounded solutions; if 7 = -|-1, Ea. (|1.3|) has quasi-periodic 
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solutions for all values of w, but we shall only consider the generic cases uj = ±1. If 70; < 0, 
in addition to plane waves and periodic waves, there exist pulse-like solutions of which are 
homoclinic as 2; ^ ibcxD to a plane wave or to the zero solution. The most famous one (if 7 = 1 
and u} = —1) is the ground state W{x) = \/2/cosh(x) which corresponds to the solitary wave of 
the focusing NLS equation. 

In contrast to the plane waves or the solitary waves which have been extensively studied [3 12H 
1^ . relatively little seems to be known about the stability of periodic waves. Spectral stability with 
respect to long- wave disturbances has been examined by Rowlands ^H]) who showed that periodic 
waves with real-valued profile are unstable in the focusing case and stable (at least in the long- 
wave regime) in the defocusing case. In the latter case, spectral stability of the small amplitude 
periodic waves has been rigorously established in using Bloch-wave analysis. Similar results 
were also obtained for certain NLS-type equations with spatially periodic potentials [HIEI- As 
for the nonlinear stability^ the only result we are aware of is due to Angulo who proved very 
recently that the family of dnoidal waves of the focusing NLS equation is orbitally stable with 
respect to perturbations which have the same period as the wave itself. We recall that the periodic 
waves of NLS with real- valued profile are called "cnoidal waves" when they have zero average over 
a period (like the Jacobian elliptic function cn), and "dnoidal waves" when they have nonzero 
average (like the elliptic function dn). 

In this paper, we study the nonlinear stability of all periodic waves of but we restrict 

ourselves to a specific class of perturbations which we now describe. Any quasi-periodic solution 
of can be written in the form 



W{x) = e'P^Qper(2A;x) , x G 



(1.4) 



where p G M, A; > 0, and Qpcr : M — > C is 27r-periodic. Here k = tt/T, where T > is the 
minimal period of \W\. The representation (|1.4|) is not unique, since we can add to p any integer 
multiple of 2k (and modify the periodic function Qper accordingly), but the Floquet multiplier e^^"^ 
is uniquely defined. Our purpose is to show that the periodic wave U{x,t) = e'(P^-^*)Qper(2A;x) 
of is stable within the class of solutions which have the same period T = n/k and the same 
Floquet multiplier e^^"^. In other words, we restrict ourselves to solutions of (jl.ip of the form 
Q^{p^-^i)Q{2kx,t), where Q{-,t) lies in the function space 



X 



i7i,,([0,27r],C) 



{q G HI^M, C) Q{z) = Q{z + 2tt) for ah z G m} 



The advantage of this restricted setting is that nonlinear stability can be established by the standard 
variational method which has been developed originally to prove the orbital stability of solitary 
waves jSl m 1^ (see also ^ Ej). However, the obvious drawback of this approach is that it 
does not give any information on the stability of the periodic waves with respect to non-periodic 
perturbations, a difficult question which remains essentially open. 

With this perspective in mind, we shall put the emphasis on the defocusing case 7 = — 1, because 
we know from ^2] that the periodic waves will be unstable in the focusing case if non-periodic 
perturbations are to be allowed. Our main result can be stated as follows: 
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Theorem 1 (Orbital stability of periodic waves in the defocusing case) 

Let 7 = —1, u) = \, and assume that W{x) = e^^^ QperC^kx) is a solution of U.cl\) with p & R, 
k > 0, and Qpcr £ ^ , «s in il-4\ )- Then there exist Co > and eo > such that, for all R ^ X 
with 1^ €0) ihe solution U{x,t) = e^^'P^~^^^Q{2kx,t) of the NLS equation Hl.l]) with initial 

data U{x,0) = e'P^(Qpcr(2A:x) + R{2kx)) satisfies, for all t E R, 

inf \\Q{;t)-e''PQpeA--mx < Co\\R\\x . (1-5) 

¥>,?e[0,27r] 

This result is known to hold for small amplitude periodic waves fI2l | , in which case the constants 
Co, eo do not depend on the wave profile Qper- Here we remove the smallness assumption, but our 
argument relies in part on the calculations made in jl2j (see Lemma 13.31 below) . 

Remarks 

1. Theorem ^includes the situation where is constant, in which case TV is a plane wave rather 
than a periodic wave. Since stability is well-known for plane waves |22j . we shall assume henceforth 
that \ W\ is a nontrivial periodic function. In such case, we emphasize that the wavenumber k which 
appears in (|1.4j) is always given by A; = vr/T, where T > is the minimal period of \W\. This is 
very important because our approach does not allow to prove the stability of periodic waves with 
respect to perturbations whose period is an integer multiple of T. 

2. It is interesting to see what Theorem ^ means in the particular case of cnoidal waves. For 
such waves we have W{x) = e^^^Qperi^kx) where p = k = tt/T and T > vr is the minimal period 
of \W\. The Floquet multiplier e^^"^ is therefore equal to —1, so that W{x + T) = —W{x), for 
all X G M. In particular, W is periodic with (minimal) period L = 2T. Theorem ^ then shows 
that the L-periodic cnoidal wave U{x, t) = e~'*T^(x) is orbitally stable with respect to L-periodic 
perturbations W provided that W{x + L/2) = —W{x) for all x S M. As explained in without 
this additional assumption the classical approach does not allow to prove the orbital stability of 
cnoidal waves with respect to perturbations which have the same period as the wave itself. 

The rest of the paper is organized as follows. In Section ^ we recall the classification of the 
bounded solutions of with 7 = — 1 and lo = 1. These solutions can be interpreted as the 
trajectories of an integrable Hamiltonian system with two degrees of freedhom, which is proved to 
be non-degenerate in the sense of KAM theory. We thereby answer a question raised by Bridges and 
Rowlands in connexion with the stability of the quasi-periodic solutions of the Ginzburg-Landau 
equation. 

Section |31 is devoted to the proof of Theorem ^ As indicated above, we follow the general 
approach of Grillakis, Shatah, and Strauss. The main difficulty is to verify the assumptions of 
the stability theorem in • We first check that the second variation of the energy functional at 
the periodic wave has exactly one negative eigenvalue. This result has been established for small 
waves in and a continuity argument allows to extend it to periodic waves of arbitrary size. 
We next consider the structure function (which is called ^^d{uj)" in [W) and show, by a direct 
calculation, that its Hessian matrix has a negative determinant. Both properties together imply 
orbital stability. 

Finally, in Section we extend our results to the focusing NLS equation. The situation is more 
complicated here, because we have families of periodic waves for all values of w G {—1; 0; 1}. By 
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comparison with the spatiahy homogeneous rotating wave U{x,t) = e'*, these periodic waves may 
be called "counter-rotating" if a; = +1, "standing" if u; = 0, and "corotating" if a; = —1. As was 
already mentioned, we restrict ourselves to the generic cases w = ±1. In both situations, we show 
that the Hamiltonian system corresponding to ()1.3() is non-degenerate in the sense of KAM, and we 
deduce as in Section 13] that the second variation of the energy functional has exactly one negative 
eigenvalue. It follows that orbital stability holds provided the Hessian matrix of the structure 
function has a negative determinant (see Propositions 14.31 and 14.61 for precise statements). Unlike 
in the defocusing case, we do not give a complete proof of this property, because the calculations are 
excessively complicated. But the determinant is easy to evaluate numerically (or even analytically 
in various parameter regimes), and it appears to be negative for all periodic waves. Assuming 
this to be true, we conclude that the analogue Theorem ^ holds in the focusing case too. Thus, 
in contrast to what happens when long- wave disturbances are considered jT^], there is apparently 
no difference between the focusing and the defocusing case as far as periodic perturbations are 
considered. 

Acknowledgements. The authors are indebted to A. De Bouard and L. Di Menza for fruitful 
discussions. This work was partially supported by the French Ministry of Research through grant 
ACI JC 1039. 



2 Properties of the periodic waves 

In this section, we study the bounded solutions of the stationary Ginzburg-Landau equation 

W^^ix) + W{x) - \W{x)\^W{x) = , (2.1) 

where : R — > C. If we interprete the spatial variable x G M as a "time", Ea. H2.1|) becomes an 
integrable Hamiltonian dynamical system with two degrees of freedhom. The conserved quantities 
are the "angular momentum" J and the "energy" E: 

J = Im(W,) , E = ^\W,\^ + ^\W\^-^\W\^ . (2.2) 

If is a solution of ()2.1|) with J ^ 0, then W{x) ^ for all x S M, so that we can introduce the 
polar coordinates W{x) = r{x) e^'^^^\ The invariants then become 

2 1^ '^r J"^ ^"^ '^^ 

J = r^io^ , E = . 

^ ' 2 2r2 2 4 

Let D C be the open region defined by 

D = |(J,^) G M2 j2 < 4/27, E^{J) <E < S+(J)} , (2.3) 

where the functions E^,Ej^ are defined in 1)2. 4() below (see also Fig. 1). Our starting point is the 
following well-known classification 131 : 
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Fig. 1: The region D in the parameter space for which Ea. l2.H has bounded solutions. 

Proposition 2.1 /f 14^ : M — > C is a bounded solution of \2.1]) . the corresponding invariants (J, E) 
belong to the closure D of D and the following properties hold: 

(i) If {J,E) E D and J ^ 0, then the modulus \W\ and the derivative of the phase of W are 
periodic with the same minimal period T{J,E) > 0. // (J,E) E D and J = 0, then W is 
real-valued, up to a phase factor, and periodic with period 2T(0,E). 

(ii) If {J,E) G dD, then either W{x) = Wg,^{x) = (1 - g2)i/2gi{gx+^) ^^^g ^ ^ 

some (p G [0,27r], or W is a homoclinic orbit connecting Wq^ip_ at x = — oo to Wq^ip^ at 
X = +00, for some q'^ <l/3 and some ip^^ip^ G [0, 27r]. 

Proof. Although the arguments are standard, we give a complete proof of Proposition 12 . II because 
it will serve as a basis for all subsequent developments in this section. Let PV^ : M ^ C be a bounded 
solution of (|2.1|) . and assume first that J / 0. Then W{x) = r{x) e^'P'^^'^ and E = ^r^ + Vj{r), 
where Vj is the "effective potential" 

^4 

^iW = ^ + y-j, r>0. 

If > 4/27, then Vj{r) < for all r > 0, hence ()2.1() has no bounded solution in that case. Thus 
we must have < 4/27. If < < 4/27, we can parametrize J in a unique way as 

J = q{l-q^) = Q{l-Q^) , where < < 1/3 < < 1 . 

Then — y + J = {y — q){y — Q){y + q + Q), hence in particular J = qQ{q + Q). With this 
parametrization, it is easy to check that Vj{r) has a (unique) local minimum at rq = (1 — Q^)^/^ 
and a (unique) local maximum at = (1 — g^)^/^ (see Fig. 2). We define 

E^{J) = y;(yr^) = i(l-Q2)(i + 3Q2) ^ (2.4) 

E+{J) = VjiVl^) = \{l-q'){l + 3q') . 

Since is a bounded solution of 1)2. II) . we necessarily have E^{J) < E < Ej^{J). This gives three 
possibilities: 
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(i) U E = E_{J), then (up to a global phase factor) W{x) = (1 — Q^)^''^ e^'^^, i.e. 14^ is a periodic 
solution with constant modulus. 

(ii) If ii^ = E^{J), then \ W\ is either constant or homoclinic to as x ^ ±00. In the first case, 
W{x) = (1 — (7^)^/^ e^"^^ (up to a phase factor). In the second case, up to a translation and a 
phase factor, we have the explicit formula 



P7(x) = (2((72+c2tanh2(cx)))'^%'''"+''^^^*""(t*'^°'^(^"» , where c = ' ' ^ 

(iii) If E_{J) < E < E+{J), the modulus r = |Vl^| and the phase derivative (fx = lm{Wx/W) 
are periodic with the same period. If we denote by ri < r2 < the three positive roots of 
E — Vj{r) as in Fig. 2, this (minimal) period is 

T{J,E) = 2 / , . (2.5) 

Jr.{J,E) ^2{E-Vj{r)) 

Another important quantity is the increment of the phase ip over a period of the modulus, 
namely 

^J,E) = 2 / 4 . = . (2.6) 

Jr,{J,E) ^2{E-Vj{r)) 

Since in general ^{J,E) is not a rational multiple of tt, the solution W{x) = r(x) e'*^^^'-* of 
(|2.1j) is not periodic, but only quasi-periodic. 




Fig. 2: When < < 4/27, the effective potential Vj has a local minimum at rq = (l — Q^)^''^ and a local 
maximum at rq = {I — g^)^''^, where q,Q are implicitely defined by J = g(l — q^) = Q{1 — Q^) and < 1/3 < . 

In the limiting case = 4/27, the effective potential Vj is stricly decreasing over M_(. with an 
inflexion point at r = (2/3)^/^. Since W is bounded, we must have E = E^{J) = E^{J) = 1/3, 
hence (up to a global phase) W{x) = (1 - q'^y/'^e^i^ with q = (l/3)^/^sign( J). 

Finally, if is a solution of (|2.1|) with J = 0, then (replacing W{x) with W{x) e"^ for some 
if E [0, 2tt]) we can assume that W{x) G M for all x G M. Then E = ^r^ + Vb(r) and we have the 
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same discussion as above with E^{0) = and -E+(0) = 1/4. U E = 0, then W = 0. If £' = 1/4, 
then either W{x) = ±1 or W{x) = it tanh(x/A/2) (up to a translation). U < E < 1/4, then W 
is periodic with half-period 

fr2{0,E) J 

r(0, ii;) = 2 / ^^^^^ , where r2(0, i?) = (1 - ^/r^i:^)!/^ . 

Jo ^/2{E - Voir)) 

Moreover, Wix + T{0, E)) = —Wix), hence \W\ is periodic with (minimal) period TiO,E). □ 

As we shall see at the end of this section, ProDosition l2.1l implies the existence of a six-parameter 
family of quasi-periodic solutions of the nonlinear Schrodinger equation (|1.1() . Before doing that, 
we study in detail the properties of the period T and the phase increment $ defined in (|2.5|) . (|2.6|) . 
because these quantities play a crucial role in the stability analysis of the solutions of (|2.1j) . both 
for the Schrodinger and the Ginzburg-Landau dynamics. We first give explicit formulas for T and 
$ (see also (HI El) which are convenient for analytical study and numerical approximation. 

Lemma 2.2 Assume that (J, E) £ D and denote by < yi < y2 < 1/3 the roots of the cubic 
polynomial Piy) = y^ — 2y^ -|- AEy — 2J^. Then 

T(J. E) = V2 /" , = 2V2 r . (") 

Jyi \/iy -yi)(y -y2)(y -y?.) Jo \/y?,-s(tp) 



where s((/3) = yi cos (99) -|- 1/2 sin ((/?). Similarly, if J ^ 0, 



mE) = V2 r i , = 2V2 r 4-. ■ 

jyi y Viy-yi)iy-y2)iy-y3) Jo ^(v?) Vys - 4^) 



1 


\ p 






y 




/yi 


yi\ 







Fig. 3: The three roots of the cubic polynomial P{y) = - 2y^ + AEy - 2.f for (J, E) e D. 

Proof. Observe that Pir"^) = Ar'^iE — Vj(r)), hence Piy) has three nonnegative roots whenever 
iJ, E) £ D (see Fig. 2). Thus, using the change of variables r = in (|2.5j) . we obtain the 
first expression in 1)2. 7|) . The last expression then follows by setting y = s ((/?), so that dy = 
2\/iy - yi)iy2 - y) dip. Similarly, if J / 0, we obtain from □ 
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It follows in particular from ()2.7|) that T is a smooth function of (J,E) £ D. In contrast $ 
cannot be extended to a continuous function over D, because (as is easily verified) 

lim <^(J, E) = ±7T , for 0<E <l/A . 
This suggests to introduce the renormalized phase ^ : D ^M. defined by 

ifJ = 0. ^ ^ 

Lemma 2.3 For any {J,E) G D, one has 

/'7r/2 J 

^iJ,E) = 2V2J / ^ ^=^^ . (2.10) 

-'0 ^/mVy3-s{lp){^/m + ^/y3-s{(p)) 

In particular, ^ : D is a smooth function. 

Proof. The integral in (|2.1U|) is a smooth function of {J,E) £ D, so it is sufficient to establish 
H2.1U() for J 7^ 0. In that case, we remark that 

r- J 1 , 2V2/TT J \ ttV2J 

2V2 / d^ = = / = 7rsign(J) , 

because 1/12/22/3 = 2J^. Thus, using (|2.8j) and (|2.9|) . we obtain 

r./2 J . I 1 



*(J,S) = 2V2 [ 

Jo 

= 2V2J 



d(^ 



d(/5 



'0 Vys Vvs - •s(v') (\/y3 + \/y3 - s((/7)) 

which is the desired formula. □ 

We next study the monotonicity properties of T and ^. Since T( J, E) is an even and ^{J, E) 
an odd function of J, we may restrict ourselves to the half-domain D-^- = {{J,E) £ D \ J > 0}. 

Proposition 2.4 

dT 

g^{J,E)>0 for all{J,E)£D. 

d"^ dT 
ii) g^iJ,E) ^ Q^iJ,E) = -QjiJ^E) > for all {J,E) G D+. 

Proof. The monotonicity of the period T with respect to the energy E has been established in 
[0]. We give here a different argument, which is also a preparation for the proof of Proposition 12. HI 
below. Let {J,E) G D. Since 2/1,2/2,2/3 are solutions of the cubic equation — 2y^ + 4:Ey — 2j'^ = 0, 
we have 2/1+2/2 + 2/3 = 2 and 

dm ^Vi dyi _ 4J --lOQ m 

dE 32/2 - 42/. + 4i? ' dJ 32/2 - 42/. + 4i? ' ^ ' ^"-''^ 
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In particular, 

^<0 ^>0 ^ + ^--^>0 (212) 
dE dE dE^ dE ~ dE ^^-^^^ 

because P'{yi) > 0, P'{y2) < 0, P'{y3) > (see Fig. 3). Similarly, if J > 0, 

^>0 ^<0 M + _M < (2 13) 

On the other hand, differentiating 1)2 .Tjl with respect to E and J, we find 

where 

l + cos2(v9) ^ ^ /- /--/^ l + sin2((^) ^ 

The crucial observation is: 

^ f-/^ sin^((^)-cos^(v.) ^ 
^2-^1 = V2/ — dip > . (2.15) 

Jo (ya - s(w)''/^ 

Indeed, this inequality follows from Lemma 12.51 below, with / = [0,7r/2], d/j, = (2/7r)d<y9, f{f) = 
sin^(v3) — cos^((^) and g{ip) = {y^ — s{(p))^'^^'^ (remark that f,g are strictly increasing, and that 
/q"/^ f{ip) dip = 0.) Thus A2> Ai> 0. We conclude that 

^ = (^2 - ^i)^ + + — j > . (2.16) 

Similarly, if J > 0, 

Finally, the relation dT/dJ = —d^/dE is a consequence of the following standard observation 
[3]: if ^ : D ^ M is the action functional defined by 



rr2(J,E) 

A{J,E) = 2 / ^/2{E-Vj{r))dr , (2.17) 



'ri(J,£) 

then a direct calculation shows that T = dA/dE and <I> = —dA/dJ, and the result follows. □ 
The following elementary result will be used several times in this paper: 

Lemma 2.5 Let n he a (Borel) probability measure on some interval / C M, and let f,g : I 
be bounded and measurable functions. Then 

I f{x)g{x)dn- U f{x)df?jU g{x)df?j = ^ / (/(a^) -/(?/)) (ff(a^) - ^(y)) d^x d/i^ . 

In particular, if both f and g are strictly increasing or strictly decreasing, and if the support of /i 
is not reduced to a single point, then 



(j{x)g{x)dfi > (^/(x)dM)(^5(x)d/x) . 
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We now establish an important non-degeneracy property of system 1)2. For any {J, E) £ D, 

let 

/ dT 9* \ 

AiJ,E) =det |J If ]iJ,E). (2.18) 
\ dJ dJ J 

Proposition 2.6 For all {J,E) e D we have A{J,E) > 0. 
Remarks 

1 . Proposition 12.61 shows that the integrable Hamiltonian system (|2.1j) is non-degenerate in the 
sense of KAM theory. Indeed, by Liouville's theorem, we can express (at least locally) the energy 
E of the system in terms of the action variables J and A, where A is defined in (|2.17|) . If we denote 
E = H{J,A), a direct calculation 5 shows that 

A - -T^det ^ ^-^^-^ 
^ - J aet g2j^ 

\ dJdA aJ^ / 
hence the KAM determinant (the Hessian of H) is always negative. 

2. By Proposition 12.41 it is clear that A > whenever d^/dJ > 0. Unfortunately, the latter 
inequality is not true for all {J,E) £ D. 

Proof. To evaluate A we have to compute the derivatives of the renormalized phase ^ with respect 
to E and J. Using (|2.1U|) . we obtain after straightforward calculations: 

dyi dvo dyi dvo 



where 



^2 Jo I yl'\y^ - s(v9))V2 (^ + v'ya - s{^)Y 

l^\y3 - siip))y^ iVm + - ^ sin2((^) J 



d<f , (2.20) 



B3 = 2V2 / dy? = - . 

(Here the formula for Bi should be read with cos^((/7) in the right-hand side, and B2 with sin^ ((/?); 
both expressions are otherwise identical.) Using 1)2. 14() . 1)2. 19|) together with the identities 



jdyi ^ dyi • 1 o Q 



which follow from 1)2. llj) . we find 



A{J,E) = {A,^-^ + A,^-^)(By-^ + B,J^-^ + B. 
^ ' ' \ dE dEJ\ dJ dJ 

a,^ + a,^)(b,j^ + b,j^) 

dJ dJ )\ dE dEJ 



dT 
^dE 



•nA.B.-A.B,{^'^-ll'^).B.. 
{A,B2 - A2B,){y, -y,)^^ + B,^ . (2.21) 
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Since dT/dE > 0, -B3 > 0, y2 — 2/1 > and {dyi/dJ){dy2/dJ) < 0, it is sufficient to prove that 

Ai =^^2^1 -^1^2 > 0. 

To achieve this goal, we set cr(y?) = (1 — and we observe that Ai = ^/2y^ '^^^Aj, 

Bi = \/2 2/3 ^^^Bj, for i = 1,2, where 

i + cos2((^) ^ ^ r/^l + sm\^) ^ 

Ai = / — :^3T^dv., A2 = / ——fld^, (2.22) 



0-((/7)3 ' cj((/9)2 

cos^(w)\l , 

1 + . 2 dV' • 

sm {if) I J 



Bi _ /"^/^j- 2 + ct((/^) ^ l + 2f7((^) 



B2 h lcT(^)(l+a(^))2 + 

Thus Ai = 2y3"'^A2, where A2 = A2B1 - A1B2. To prove that A2 > 0, we remark that 

1 + 2(7 1 1 ,2 + (7 l + cos2(o9) o + ^\Y?(ip) 

— and — 



a3(l + f7)2 a3 a(l + a)2 ' a(l + a)2 a{l + a)^ a{l + ' 

It follows that Bi = Ai + Bi and B2 = A2 + B2, where 

- _ aM+sin^M ~ _ a{^) + co.\^) 

Now, we know from H2.15() that A2 > Ai > 0, and the same argument (using Lemma l2 .51) shows 
that Bi > B2 > 0. Thus 

A2 = A2B1 - A1B2 = A2B1 - A1B2 > , 
and the proof is complete. □ 




Propositions 12 .41 and 12 .61 implv that the quasi-periodic solutions of 1)2. 1|) can be parametrized by 
the period T and the renormalized phase ^, instead of the angular momentum J and the energy 
E. Indeed, let 



D 



^{T,^) T > TT , \^\ <^{T)^ , where ^(T) 



T2 + 27rSi/2 

— vr . 
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Proposition 2.7 The map (T, : D ^ D is a smooth diffeomorphism. 

Proof. Let J = Q{1 — Q^), where 1/3 < < 1. Using (|2.7|) it is straightforward to verify that 
lim T(J, E) = = , , lim Tf J, E) = +oo . (2.24) 

Similarly, if l/\/3 < Q < 1 (so that J > 0), we deduce from (|2.<S|) that 

lim $(J,S) = ^^'^ , lim $(J,^) = +00 . (2.25) 
E^E.{.j) V3g2 _ 1 i?^£;+(J) 

Since dT/dE > by Proposition |231 it follows from (|2.24|) that the range of the map T : D ^ M is 
exactly the interval (vr, +oo). Fix Tq > vr and let S = {(J, E) & D \ T{J, E) = Tq}. By the Implicit 
Function Theorem, S is a smooth curve in D which can be represented as a graph over the J-axis. 
Moreover, we know that S connects the boundary points {—Jq, E^{Jq)) and (Jq, -E'_(Jo)), where 
Jo = Qo{^ — Qo) ^-iid Qo G (l/\/3, 1) is determined by the relation 

7r\/2 _ ^ 

(see (|2.24p ). Now, Proposition 12.61 implies that the restriction of ^ to the curve E is a strictly 
increasing function of J, because 

dJ ~ KdEJ KdE'dJ ~ dJdEJ ^ ' 
Thus ^ varies from — ^'o to ^'o on the curve S, where by ()2.25|) 

vrV2Qo /r2 + 2vr2\i/2 



This proves that (T, ^) : — > D is onto, and the monotonicity properties established in Proposi- 
tion E31i™ply that (T, ^') is also one-to-one. □ 

To conclude this section, we briefly verify that Proposition 12.11 implies the existence of periodic 
waves of the NLS equation of the form (|1.2|) . Fix {J,E) E D, and let : M — > C be a bounded 
solution of ()2.1|) satisfying (|2.2|) . (The proof of Proposition 12.11 shows that this solution is unique 
up to a translation and a phase factor. If needed, we can specify a particular solution by imposing 
for instance VF(0) = r2, VF'(O) = iJ/r2, where r2 > is as in Fig. 2.) We now set 

Wix) = e'^^ P{kx) , a; G R , (2.26) 

where 

k = , and £ = Zii'^l . (2.27) 

T{J,E) T{J,E) ^ ' 

As X I— > |VK(x)| is periodic with minimal period T{J, E), it is clear that y \P{y)\ is periodic with 
minimal period vr. Moreover, since W{x + T) = e'^W{x) by definition of <^{J,E), we also have 
P[y + vr) = -P{y) for all y e M, hence P is 27r-periodic. Thus U{x, t) = e-'^W{x) = e'(^^-*)P(/cx) 
is a quasi-periodic solution of (|1.1|) of the form (|1.2|) . with uj = 1 and c = 0. 
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Remark 2.8 Using the continuous symmetries of the NLS equation, we can produce for each pair 
{J,E) G D a four-parameter family of periodic waves: 

Ucx^,^{x,t) = \SP^'-^-^^^^'-^^P{k\{x-ct)+i) , 

where c G M, A > 0, (/?, ^ G [0, 27r], and p^^x = X£ + c/2, oj^^x = + cXi + Taking into account 

the parameters J,E, we obtain altogether a six-parameter family of periodic waves of 

The representation l|2.26jl is well-adapted to understand the connection between the bounded 
solutions W of H2.1() and the periodic waves of the NLS equation, especially in the case of the 
cnoidal waves for which J = i = 0. However it is more convenient for our purposes to write the 
solution W of H2.ll) in the alternative form 

W{x) = e'^^+''>Q+{2kx) = e'(^-'=)^Q^(2A;x) , x G M , (2.28) 

where Q^{z) = e^^^^'^P{z/2). By construction, and \Q^\ are now periodic functions with 
the same minimal period 27r. This property facilitates the description of the special class of 
perturbations that we use for the stability analysis, see the statement of Theorem ^ and the 
remarks thereafter. The representation ()2.28() is also very natural for solutions which are close to 
plane waves: in such a situation, either or Q~ is close to a constant, depending on the sign of 
J. This follows from the fact that ^ + A; = ^/T if J > and £ - /c = $/r if J < 0. 



3 Orbital stability 

Our aim in this section is to show that the periodic waves of the defocusing nonlinear Schrodinger 
equation Hl.l|) with 7 = — 1 are stable within the class of solutions which have the same period T 
and the same Floquet multiplier e'*. Given {J,E) G D, where D C is the parameter domain 
(|2.3|) . we consider the periodic wave Uj^E{x,t) = e~^^Wj^E{x), where Wj^e is the unique solution 
of (|2.1|) with initial data Wj^e{0) = ?'2, Wj ^{0) = iJ/r2, and r2 > is as in Fig. 2. In particular, 
Wj^E satisfies Let p = k + £, where k = kj^E and £ = £j^e are defined by ((T77|) . As in (gUHl), 

we set Wj^Eix) = e^^^Qj^Ei'^kx) and we recall that Qj^e{z) is a 27r-periodic function of z. 

To study the stability of the periodic wave U j^e{x, t), we consider solutions of (|1.1|) of the form 

U{x,t) = e'(?'^-*)Q(2fcx,i) , (3.1) 

where Q{z,t) satisfies the evolution equation 

iQt + 4ipkQ, + ik'^Q,, + (1 - /)Q - |Q|2q = . (3.2) 

By construction, Qj^e{z) is now a stationary solution of 1)3.2(1 and our goal is to show that this 
equilibrium is stable with respect to 27r-periodic perturbations. We thus introduce the function 
space 

X = i?j^g,([0,27r],C) = |n G Fic(M,C) u{z) = u{z + 2tt) for all z G m} , 
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which is viewed as a real Hilbert space equipped with the scalar product 

f-2TT 

{u,v)x = Re / {u{z)v{z) + Uz{z)vz{z)) dz , u,v G X . 
Jo 

As usual, the dual space X* will be identified with Hpj^{[0, 27r], C) through the pairing 

f-2lT 

{u,v) = Re / u{z)v{z) dz , u G X* , v X . 
Jo 

It is well-known that the Cauchy problem for is globally well-posed in the space X (see 
[51 113L HHl I17j ). Moreover, the evolution defined by (|3.2|1 on X is invariant under the action of a 
two- parameter group of isometries: the space translations and the phase rotations. The symmetry 
group G is thus the two-dimensional torus T-^ = (]R/27rZ)^ which acts on X through the unitary 
representation TZ defined by 

(7^(^,5)n)(z) = e-'^^uiz + , ueX , {ip,O^G. 

In view of these symmetries, it is natural to introduce the semi-distance p on X defined by 

p(u,v) = inf \\u — Tl(,„ f:\v\\x , u,v G X . (3.3) 

The main result of this section is the following reformulation of Theorem ^ 

Proposition 3.1 Given {J,E) G D, there exist Cq > and eo > such that, if Qq £ X satisfies 
p{Qo,Qj,e) < e for some e < eo, then the solution Q{z,t) of iV. i^) with initial data Qq satisfies 
p{Q{;t),Qj,E)<Coefor alltGR. 

To prove Proposition 13.11 we follow the general approach of Grillakis, Shatah, and Strauss 
jl5l I16j . We first observe that Eq. (|3.2|) inherits from the original NLS equation several conserved 
quantities: 



1 



27r 



r2TT 



M{Q) ='-jJ Q{z)Qz{z)dz , (3.4) 

£{Q) = fj (2k'\QM\'' + \\Q{zt)<^z , 

which will be referred to as the charge, the momentum, and the energy, respectively. The charge 
N is conserved due to the phase invariance of 1)3. 2(1 . the momentum M due to the translation 
invariance, and the energy £ because ()3.2() is autonomous. Clearly N , M, and £ are smooth real 
functions on X. Their first order derivatives are therefore smooth maps from X into X*: 

N'{Q) = Q , M'iQ) = \Qz , £'{Q) = -Ak^Q,, + IQpQ . (3.5) 
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Similarly, the second order derivatives are smooth maps from X into C{X,X*), the space of all 
bounded linear operators from X into X*: 

N"{Q) = 1 , M"{Q) = Id, , 8"{Q) = -4k^d,, + IQI^ + 2Q ® Q , 

where {{Q Q)u, v) = J^^ Ke{Qu) Ke{Qv) dz for all u,v G X. 

By construction, the periodic wave profile Qj^e is a critical point of the modified energy 

SjMQ) = £iQ) - (1 - P^)N{Q) - ApkM{Q) , (3.6) 

namely £j ^{Qj^e) = 0. To determine the nature of this critical point, we consider the second 
variation 

Hj,E = £j,EiQj,E) = -^k^d,, - Aipkd, - (1 V) + \Qj,e\^ + 2Qj,E Qj,E ■ 

Since X is a real Hilbert space, it is natural to decompose its elements (which are complex functions) 
into real and imaginary parts, in which case we obtain the matrix operator 

^ ^ / -4^8,, - (l-p2) + 3RIe + IIe ^Pkdz + 2Rj,eIj,e 

■^'^ V -Apkd, + 2Rj^Elj,E ' -Ak^d,,-{l-p') + RlE + ^ll 

where Qj^e = R,j,e + iI,j,E- As is easily verified, Hj^e is a self-adjoint operator in Xq := 
Lpgj,([0, 27r], C) with compact resolvent, and Hj^e is bounded from below. According to ^lEl, a 
crucial information is the number of negative eigenvalues of this operator. 

Proposition 3.2 For any (J, E) € D, the operator Hj^e acting on Xq has a simple negative 
eigenvalue, a double eigenvalue at zero, and the rest of the spectrum is strictly positive. 

Proof. Since Ea. ()3.2() is invariant under the action of G, we know that TZ(^^^^^Qj^e is a stationary 
solution for all {(f,^) S G. Differentiating with respect to and ^ at (</?,^) = (0,0), we obtain 
Hj^eQ'j E — H,j^e{}Q,j,e) = 0. It is clear that the functions Q'j ^ and '\Qj^e are linearly independent 
in the domain X2 := Hp^j.{[0, 27r], C) of Hj^e (otherwise Wj^e{x) = e^^^Qj^EC^kx) would be a plane 
wave, contradicting the assumption {J,E) G D). Thus we see that zero is an eigenvalue of Hj^e of 
multiplicity at least two. 

Following [9, , we next show that the multiplicity of zero as an eigenvalue of Hj e in Xq is always 
exactly two. The idea is to produce two other solutions of the differential equation HjeQ = by 
differentiating the profile Qj^e with respect to J and and to verify that none of these solutions 
belongs to the domain X2. Since the kernel K = {Q G C'^(M, C) | Hj^eQ = 0} is a four-dimensional 
(real) vector space, it will follow that Kr\X2 = spanjQj^, i(5j,_E}, which is the desired result. 

However, the profile Qj^e is a stationary solution of the equation (|3.2() which has coefficients 
depending upon J and E, so that we cannot find the two solutions of HjeQ = just by dif- 
ferentiating Qj^E with respect to J and E. Instead, we start with the quasi-periodic solution 
Wj^e{x) = e^^^Qj^E^^kx) of (|2.1|) . In view of Proposition 12. 7[ we can use T and ^ (instead of J 
and E) to parametrize the family of quasi-periodic solutions of (|2.1j) . For our present purposes, 
the most convenient set of parameters will be {k,p), where p = k + £ and k,i are defined in (|2.27|1 . 
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Since the equation l\2.1\i does not depend upon k and p, the derivatives dWj^E/dk and dWj^E/dp 
are solutions of the hnear equation 



io" + UJ- 2\Wj,E\'^io - WIe^ = . 
It is then straightforward to check that the functions Ri{z) and R2{z) defined for z = 2kx through 

Ri{2kx) = e-'v-^^^[x) = ^^{2kx) + 2xQ'j^E{2kx) , (3.7) 

R2{2kx) = e-'P-^^^{x) = ^^{2kx)+ixQj^E{2kx) , 
op op 

satisfy HjeRi = Hj^eR2 = 0, namely Ri,R2 G K. Of course -Ri,i?2 depend on the point 
(J, E) £ D where the derivative is taken, but we omit this dependence for notational simplicity. 

Now, let R{z) = AiRi{z) + A2R2iz) for some ^1,^2 G ^- Using the definitions and the 
fact that Qj^E is 27r-periodic we find that R satisfies the periodicity conditions R{2tt) = R{0) and 
R'{2Tr) = R'{0) if and only if 

2AiQ'j^EiO) + -iA2Qj,Em = , and 2^iQ'j_^(0) + iAaQ'j.sW = . 
This linear system has a nontrivial solution (^1,742) G if and only if 

Q'M^f = QjMo)Qj,Eio) . (3.8) 



Using again the quasi-periodic solution Wj^Eix) which satisfies the initial conditions Wj^e{0) = ^2, 
Wje{0) = \J/r2 and WJ e{^) = — r2, it is not difficult to verify that 1)3. 8|) is equivalent to 
j2 = rl — r^. But we know that J"^ = fq — t^q = — with rq < r2 < Vq (see Fig. 2), hence 
J2 < r| - so that (ESJ never holds. We conclude that R = AiRi + A2R2 belongs to X2 only 
if Ai = A2 = 0. This shows that the four functions Q'j eAQj,e, Ri, R2 are linearly independent 
(over M) and therefore form a basis of the kernel K. Moreover, K Ci X2 = spanjQj^, iQj^^;} as 
expected. 

It is now easy to conclude the proof of Proposition 13.21 Indeed, since the eigenvalues of Hj^e 
depend continuously on {J,E) G D, the fact that the zero eigenvalue has constant multiplicity 
implies that the number of negative eigenvalues of Hj^e remains unchanged when (J, E) varies 
over D. It is therefore sufficient to verify that the conclusion of Proposition 13.21 holds for one 
value of (J, E) G D. But as (J, i?) — > the wave profile Qj,_e converges uniformly to zero, hence 
Hj^E converges to a constant coefficient operator whose spectrum is easily determined by Fourier 
analysis. Using a perturbation argument, it is then straightforward to localize the eigenvalues of 
Hj^E for (JjE) close to the origin. This calculation is performed in jl2j . using an appropriate 
parametrization of the small amplitude periodic waves. In the particular case J = 0, we obtain 
the result below, which concludes the proof of Proposition 13.21 □ 

Lemma 3.3 f | 12l Proposition A.l]) If J = and E > is sufficiently small, the operator Hj^e 
acting on Xq has exactly three eigenvalues {Ai,A2,0} in a neighborhood of the origin, and the rest 
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of its spectrum is positive and bounded away from zero. Moreover, the eigenvalues Ai, A2 are simple 
and satisfy ^ 

\i{E) = -E + 0{E'^) , X2{E) = 3E + 0{E'^) , as E ^ . 

Remark 3.4 By construction the wave profile M^j,_e(x) = e^^^Qj^E^'^kx) is a stationary solution 
of the time- dependent Ginzhurg- Landau equation 

Wt{x,t) = W,,{x,t) + W{x,t)-\W{x,t)\''W{x,t) . (3.9) 

If we look for solutions of h'J. y\) of the form W{x,t) = e^P^Q{2kx,t), we obtain the evolution equation 

Qt{z, t) = ik^QU^, t) + 4ipkQ,{z, t) + (1 - p^)Qiz, t) - \Q{z, t)pQ(z, t) , (3.10) 

which is very similar to In particular, the linearization of \'3.1U\) at the equilibrium Qj^e 

is Qt = —Hj eQ- Thus Provosition VJ.'A imvlies that the quasi-periodic solutions of the Ginzburg- 
Landau equation y\) are always unstable, even within the class of solutions with the same period 
and the same Floquet multiplier as the original wave profile. This complements previous results by 
Bridges & Rowlands /5/, and by Doelman, Gardner & Jones which show that the quasi-periodic 
solutions of \S. are unstable with respect to long-wave disturbances ("sideband instability" ) . 

We now continue with the proof of Proposition 13.11 The next observation is that, for any 
(J, E) S D, the equilibrium Qj^e of (|3.^B is a member of a two-parameter family of travelling and 
rotating waves of the form 

Q{z, t) = e-'""* Q^E^^ + ct) , z G M , t G M , (3.11) 

where {uj,c) lies in a neighborhood of the origin in (the Lie algebra of G), and the profile 
Qjg G X is a smooth function of {oj, c) with Qj^ = Qj,E- (Actually we even have a four-parameter 
family of such waves if we take into account the action of the symmetry group G.) Indeed, take 
{J',E') e D close to {J,E), and define k', f by the formulas (I^TTj) with (J, E) replaced by {J',E'). 
Then 

U{x, t) = e^(p'^-*) Qj,^E'{2k'x) , where p' = k' + £' , 

is a solution of (|1.H) . but it is not of the form (|3.H) because p' ^ p and k' ^ k in general. However 
we can transform U{x,t) into a solution of (jLlj) of the form ()3.1() . 1)3. 11() by applying successively 
a dilation of factor A and a Galilean transformation of speed v, where 

A = ^, v = 2{Xp'-p). (3.12) 

After some elementary algebra, we obtain Q^'Ei^) = XQj'^E'i^) with 

uj = X\l - p'^) - {1 - p^) , c = AX^k'p' -4kp . (3.13) 

^The expansions of Ai, A2 do not appear explicitly in |12| . but they follow easily from the proof of Proposition A.l 
there if one observes that the case J — corresponds to a = 6 and E = 2a^ + O(a^) in the notation of |12|. 
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The mapping {J',E') (w,c) that we have just defined is a diffeomorphism from a neighbor- 
hood of {J, E) onto a neighborhood of (0,0). Indeed, using (|2.27|) and ()3.12|1 . the formulas H3.13|) 
can be written more exphcitly as 



UJ 



V2-(^ + ^')') > C = —{TT + ^')-4kp , (3.14) 



where T' = T(J',E') and ^' = ^{J',E'). Differentiating these expressions with respect to J',E' 
we obtain 



du! dc 



duj dc 
dJ' dJ' 



\A 4§I / dE' dE' \ _ TT^dE l^dE' TV dE \ /■? i c:\ 

{J',E')={J,E) ^—(dJ-Pdj) —dJ 



By Proposition 12.61 we have det(A^j^£;) = {Sk"^ /■k'^)A{J, E) > 0, hence the matrix Ai j^e is invert- 
ible. This proves the existence of the travehing and rotating waves (|3.11|) for sufficiently small 

(w,c) e M^. 

By construction, the profile Q'j''^ G X is a critical point of the functional 

£j:e(Q) = £jMQ) - ^N{Q) - cM{Q) , QeX . (3.16) 
For later use, we remark that, if A is as in (|3.12)) . we have the identity 

Sj'^iXQ) = X'^Sj'^E'iQ) , for all Q G X . (3.17) 
Following we now define, for (cj,c) in a neighborhood of (0,0), 

dj,E{u^,c) = £j;EiQTE) ^ >^^Sj',E'iQj',E') ■ (3.18) 

A crucial role in the orbital stability argument will be played by the Hessian matrix of dj e, namely: 



E 



~duj^ duj dc 
dcdLu dc? 



((^,c)={0,0) 

Proposition 3.5 For all {J,E) e D we have det(TCj^E) < 0. 

Remark 3.6 Since TCj^e is a symmetric matrix, Provosition 15'. 51 imvlies of course that TCj^e has 
one positive and one negative eigenvalue. This can be verified explicitly in the case of small ampli- 
tude periodic waves, for which the following asymptotic result is established in 

^J,i? = |(l^ ~^^{1 + 0{E)), as {J,E) ^{0,0) inD . 

Since the proof of Proposition 13 . 51 is long and technical, we postpone it to the end of this section, 
and we now show how Propositions 13.21 and 13.51 together imply Proposition 13.11 The arguments 
here are rather classical and can be found in 16 , so we shall just indicate how the general theory 
of applies to present case. 

For any (J, E) £ D, Proposition 13.21 shows that the wave profile Qj,_e is a degenerate saddle 
point of the energy Sj^e, with one unstable and two neutral directions (the latter are due to the 
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fact that (|3.2() is G-invariant). To get rid of the unstable direction we remark that the evolution 
of ()3.2|) does not take place in the whole function space X, but on codimension two surfaces where 
the charge and the momentum M are constant. Let 



^J,E = {Q£x\ NiQ) = N{Qj^e) , M{Q) = M{Qj^e)] ■ 



(3.19) 



It is easy to verify that Tij^e is indeed a smooth submanifold of X of codimension two in a 
neighborhood of Qj^e-, or more generally in a neighborhood of the orbit of Qj^e under G. The 
crucial point is that the functional £j^e is coercive on with respect to the semi-distance p 
defined by 1)3. 3|) . More precisely, there exist positive constants Ci,5 (depending on J^E) such that 

£jAQ)-£j,E{Q.hE) > Cip{Q,Q.j,E? , 



(3.20) 



for all Q e Tjj^e such that p{Q, Q,j,e) < 5- 

To prove (|3.20j) we first note that the tangent space to Sj,_e; at Qj^e is 



Tj,E = {q G X I Q) = {iQ'j^E, Q)=0} , 



(see H3.5() ). We next introduce the "normal space" Mj^e = spanlduiQj^EjdcQj^E}, where 



d 



J,E 



(a.,c)=(0,0) 



dcQj,E = ttQ^j'^ 
oc 



J,E 



((^,c)={0,0) 



As Qj'e is a critical point of £jei we have £'j e^Qj'e) ~ ^^'{Qj'e) + '^^'(Qje)- Differentiating 
this relation with respect to uj and c at (w,c) = (0,0), we obtain 

Hj,E{du.Qj,E) = N'iQj^ = Qj^E , HjAd,Qj,E) = M'iQj^ = iQ',^ , 



hence {Hj^eQi,Q2) 
find 



for all Qi G J^j,E, Q2 £ Z/.s- In a similar way, differentiating (|3.1(S|1 we 



J,EI ' 



(3.21) 



hence 



T~(-J,E 



ImTE) 



dc 



J, El 



WcM{QTe) 



(3.22) 
(3.23) 



_ _ f {Hj,E{d^Qj,E),duQj,E) {Hj,EidcQj,E),d^Qj^E) 
V {Hj^E{d^Q.J,E),dcQj,E) {Hj,E{dcQj,E),dcQj,E) 

Thus —Hj^E is the matrix of the bilinear form bj^EiQi, Q2) = {Hj^eQi, Q2) restricted to the sub- 
space Mj^E- As TCj^E is non-degenerate by Proposition 13.51 we have Mj^e H Tje = {0}, hence also 
X = Mj^E © Tj^E- Moreover, since Hj^e has only one negative eigenvalue by Proposition 13. 2| and 
since the restriction of hj^E to Mj^e has a negative eigenvalue by Proposition 13.51 a standard signa- 
ture argument shows that the restriction of hj^E to the complementary space Tj^e is nonnegative. 
More precisely, if we decompose Tj^e = Kj^e © "^j e where Kj^e = spanjQj^;, ^Q,j,e} is the kernel 
of Hj^E in X and Tj ^ = {Q G Tj^e \ {Q'j e->Q) ~ {^Qj,E,Q) = 0}, then there exists C2 > such 
that 

{Hj,eQ, Q) > C2\\Q\\l , for all Q G Tj^e ■ (3-24) 
Estimate (j3.20p follows from its infinitesimal version (|3.24j) using a Taylor expansion of the energy 
£j^E at Qj^E and the fact that £j^e is G-invariant (see JH] for details). 
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Remark 3.7 The argument above does not give any lower bound on the constants Ci, C2 in h'j.2U\) . 
\3.24\j - In particular, we cannot take for C2 the lowest positive eigenvalue of the operator Hj^e, 
unless the subspace Nj^e contains precisely the eigenfunction associated to the (unique) negative 
eigenvalue of Hj^e- However, one can check that the constants Ci,6 in \S.2i^) are hounded away 
from zero uniformly for {J,E) in any compact subdomain of D. 

With inequality H3.2U() at hand, it is now easy to conclude the proof of Proposition |23 Assume 
that Qq £ X satisfies p{Qo,Qj,e) < e for some small e > 0. Replacing Qq by TZ^^^^^^Qq if needed, 
we can assume that \\Qo — Qj^eWx < £• As Qj^e is a critical point of Sj^e, we have \£j,EiQo) — 
£.j,e{Qj,e)\ < C'se^, for some C3 > (depending on J,E). We distinguish two cases: 

1. If Qq G Tjj^Ei then the solution of (|3.2|) with initial data Qq satisfies Q{-,t) E for all t G R, 
and (|3.2())) implies that 

Cip{Q{;t),Qj,Ef < £j,E{Q{-,t))-8j,E{Qj,E) = £j,E{Qa)-£j,E{Qj,E) < Cse^ , 

for all t G M, provided e is small enough so that C^e^ < CiJ^. This is the desired result. 

2. The case where Qq ^ ^j,e can be reduced to the previous one by the following argument. 
In view of ()3.22|) and Proposition 13.51 the map (<jJ,c) ^ {N (Q^j'j^) , M {Q^j''^)) is a diffeomorphism 
from a neighborhood of (0,0) onto a neighborhood of {N{Qj^e),M{Qj^e))- Thus there exists a 
unique (w,c) G with \lo\ + |c| = 0{e) such that N{Qq) = A^(Qj;^), M{Qq) = MiQ'^^f^). By 
construction Q^/e = ^Qj',E', where {J',E') is e-close to {J,E) and A is as in 1)3. 12() . Thus, if we 
denote Q{z,t) = \~'^Q{z,t), we see that Q{-,t) G S j/^e' for all t G M. Using ()3.17|) together with 
(|3.2U() . we thus obtain 

£j:EiQi-^t))-£j:E(Q7E) = >^'{£j',E'{Q{;t))-£jMQj',E')) 

> C[\^p{Q{;t),Qj,E'f = C[X^p{Q{-,t),Q%f , 
as long as p{Q{-,t),Qy^) stays sufficiently small. Since 

£j:E{Q{-^t))-Sj:Em,E) = ^tEm - ^Te^qTe) < c'sWQo - qTeWx , 

and IIQo - Q'j'^eWx < \\Qo - Qj,e\\x + \\Qj,e - Qj^slU < C'4e, we conclude that 

p(.Q{;t),Qj,E) < p{Q{-,t),Qyi) + \\Q% - Qj^eWx < C^e for alH G M , 

if e > is sufficiently small. This concludes the proof of Proposition 13.11 □ 

Proof of Proposition 13.51 We start from the expression (|3.22|) of Hj^e- We recall that {to, c) can 
be parametrized (in a neighborhood of the origin) by {J',E') according to the formulas 1)3. 13() or 
1)3. 14() . and that Q^e ~ ^Qj',E' where A is given by ()3.12|) . Thus Hj^e = —{M j,e)~^K,j,e, where 
Mj^E is defined in 1)3. 15() and 



WE WW\ def / m'N{\Qji^E') ■^M{\Qj>^e' 



8'N d'M / \ d ATr\r^ \ 9 



d'J d'J ^ \ QJ'^" V'^J',E' J dJ. 



tN{XQj>e') 7fFM{XQj',E' 



{J',E')={J,E) 



Since we already know that det{M.j^E) > 0, it remains to verify that det(/Cj^£;) < 0. 
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We first give more explicit formulas for the coefficients of ICj^e- As Wj^e{x) = e^^^Qj^EC^kx), 
where p = k + £ and k, i are given by H2.27p . it follows from H3.4|) that 

N{Qj,e) = I r \Qj,E{z)\''dz = k r \Wj,e{x)\'' dx , 
^ Jo Jo 

M{Qj,e) = \j^Qj,E{z)Q'jAz)dz 

= \ {Wj,E{x)W'j^E{x)-ip\Wj,E{x)\'')dx = ^N{Qj^e)-i^JT . (3.25) 



Moreover, proceeding as in Lemma 12.21 we find 

N{Qj,e) = 2k / dr = 2V2k ^"^^ dip , (3.26) 

Jn{j,E) y/2{E-Vj{r)) Jo V^s - s(v') 

where s{ip) = yicos^((/9) + y2sm'^{ip) and yi,y2,y3 are as in Fig. 3. Of course, since N,M are 
quadratic functionals, we have 

N{XQj',E') = X^NiQj,E') = ^NiQj,E') , M{XQj,E') = ^M{Qj,e') ■ (3.27) 

Differentiating (|3.27jl with respect to J', E' and using (|3.25|) . (|3.26|) . we obtain after straightforward 
calculations: 

We - ^'dE^^'dE^^'dE' 5^7 - ^^a/ +^^97 + ^^97 ' ^^'^^^ 

d'M _ p d'N N 3 J dT d'M _ p d'N N 5^- 3J dT T 

'We ~ 2kWE ^ l^dE ~ YdE ' 'WJ ~ 2fc 5^7 2^ 97 ~ Ta7 ~ 2" ' 

where = N{Qj^e) and 

Cs = 2^2- ^M=d^ = ^N{Qj,e) . 

We now compute det(/Cj^£;). Replacing the second column of the matrix K,j^e by its expression 
(|3.29|) . we see that det(/Cj^£;) = Ai + A2 + A3, where 

'NB-i T\d'N 
Ai ' ^ ^ 



2tt 2 J d'E ' 

^^'^ AT / <M. 9^ \ n T / d'N dT 

Here ^3 is as in ^TW . 
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We claim that Ai < 0. Indeed, arguing as in (|2.15|1 and using Lemma 12.51 we find 

C2-Ci = V2k '^^^ ~ (sin2(c^) - cos2(vp)) d(^ > , 

because ip ^ {2y^ — s{(p)){y^ — 5(92))"^/^ is strictly increasing over [0, 7r/2]. Thus C2 > Ci > 0, 
and proceeding as in (|2.16j) we deduce that Ci{dyi/dE) + C2{dy2/dE) > 0. Since C3 > and 
dT/dE > by Proposition l2.41 we conclude that d'N/d'E > 0. On the other hand, since s{ip) < y2 
for all if G [0,27r], it follows from (|T7I), (|T77|l . (IT^ that iV < 7ry2- Moreover, as 



2/2 m ViE - Vvs - y2 

— , — / , — , / \ J — ^ ' 



it follows from ^1^, l^TM that B3y2 < T. Thus NB^ < vrT, hence Ai < 0. 
It remains to verify that A2 + A3 < 0. By ((TTH) . we have 

and similarly for the J-derivative. Thus using H2.19() and proceeding as in ()2.21|1 we find 
A2 = ^{y2-yi)jjjj{{Ci + CsA,)B2-{C2 + C3A2)B,) . 

In a similar way 

As = -|{y2-y,)'^^{C^A2-C2A,). 

Thus A2 + A3 = (2/2 - yi){dyi/dJ){dy2/dJ)A4, where 

o MC N 

A4 = -iAiC2 - A2C1) + -^iAiB2 - A2B1) + —iCiB2 - C2B1) . 

As y2 > yi and {dyi/ dJ){dy2/ dJ) < 0, we have to verify that A4 > 0. 

Like in the proof of Proposition 12.61 it is convenient to express the various constants Ai,Bi, Ci 
in terms of the new function cr{ip) = (1 — y'^^s{ip))^^'^ . Then Ai = V^y^ ^^^Aj, Bi = V^y^ ^''^Bj, 
Ci = V2ky^ ^^^Cj, where Aj, Bj are defined in 1)2.22(1 and 



2cos2((^) + (l-a( v?)^)sin^M 



Cl = / I7TN3 ^(^ ' 



2sin^(9.) + (l-aM^)cos^M 
= I 



Thus 



A4 = ^ f l(AiC2 - A2C1) + -^(AiB2 - A2B1) + -^(CiB2 - C2B 

Unfortunately, the three terms in the right-hand side are not all positive. Indeed, we recall that 
Bj = Aj + Bj for i = 1, 2, where Bi, B2 are defined in (|2.23() . Similarly, Cj = Aj — Cj, where 



-/2sin2(^ ~ r/'cofi^ 

Ci = / — -— d(p , C2 = — dip . (3.30) 
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We also recall that A2 > Ai, B2 > Bi, C2 > Ci, but Bi > B2, Ci > C2. All these bounds follow 
from Lemma I2. 51 as in ()2.15|) . Thus 

A1C2-A2C1 = CiA2-C2Ai>0, 

A1B2 - A2B1 = A1B2 - A2B1 < , (3.31) 
C1B2 - C2B1 = (A1B2 - A2B1) + (C2B1 - C1B2) < . 

As TV < 7ry2 < t^Vz and C3 = {kN/i:) < kys, it follows from (|3.31|) that A4 > (2fc/y|)A5 where 
As = ^(AiC2 - A2C1) + ^(AiB2 - A2B1) + ^(CiB2 - C2B1) . 

It remains to verify that A5 > 0. Since Bj = Aj + Bj and Cj = Aj — Cj for i = 1,2, we have 
As = (CiA2-C2Ai) + (AiB2- A2Bi) + i(BiC2-B2Ci) 

A2 - Ai + ^'~^' ) (Ci - Bi) + (Ai - ^) ((Ci - C2) - (Bi - B2; 

We claim that all terms in the right-hand side are now positive. Indeed, we already know that 
A2 - Ai > and Ci - 62 > 0. Using (IT^ and (HUpi we obtain 

~ ~ (2 + aM)sin^M-l ^ sin^((^)-cos^(v.) ^ 

because 99 1— > (1 + a{ip))~'^ is strictly increasing on [0,7r/2]. On the other hand, since (j{(p) < 1 we 
have 



^1 > / > / — = Ci , 



hence Ai — Ci/2 > a fortiori. Finally, 

2 + ^(99) 



(Ci - C2) - (Bi - B2) = r ^'^^''l^^ (sin2(v.) - cos2((^)) d(/. > , 

because 1-^ (2 + (t((/3))/(1 + cr{ip))'^ is strictly increasing on [0, 7r/2]. Thus As > 0, and the proof 
of Proposition 13.51 fhence of Proposition 13.11) is now complete. □ 



4 The focusing NLS case 

In this final section we show how the preceding results can be extended to the focusing NLS 
equation (jl.lj) with 7=1. In this case, Ea. H1.3() has quasi-periodic solutions for all values of w, 
but we restrict ourselves to the generic cases lo = \ and a; = — 1, which we consider separately. 
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4.1 Counter-rotating waves {lu = 1) 

Proceeding as in Section [21 we first study the solutions of the stationary equation 

W^^{x) + W{x) + \W{x)\'^W{x) = , x£R. (4.1) 

This is again an integrable Hamiltonian system with conserved quantities 

J = Im(W,.) , E = ^\W,\^ + ^\W\^ + j\W\^ . (4.2) 

In particular, the effective potential Vj(r) = J^/(2r^) + + r^/4 is now strictly convex for any 
J G M. If we set J = q{q^ — 1) where g G M, \q\ > 1, then the unique minimum of Vj is attained 
at r = Tq = q'^ — I and has the value 

E^{J) = Vj{y^^) = \{q' - l){3q' + 1) . 
It follows that 1)4. 1(1 has quasi-periodic solutions if and only ii {J, E) £ D, where 

D = \^{J,E) eR^ E> S-(J)} . 
The period T and the phase increment ^ of these solutions are given by 

T( J, E) = 2V2 r'' -r^= , HJ, E) = 2V2 -r== ' (4-3) 

where s{if) = yi cos^(c/?) + y2sin^((/?) and < < yi < y2 are the roots of the cubic polynomial 
P{y) = —y^ — 2y^ + 4£'y — 2J^. Similarly, the renormalized phase (|2.9|1 satisfies 

^{J,E) = -2V2J '^'^ ^ . (4.4) 

Jo yf^\/s{^) -y-i{^f^+ ^Js{ip) -yz) 

In contrast with the defocusing case, the period T( J, E) is now a decreasing function of the 
energy. The analogue of Proposition [TH is: 

Proposition 4.1 

dT 

^) g^{J,E)<0 forall{J,E)eD. 

d"^ dT 
ii) -Q^iJ^E) = —QjiJ,E) > for all {J,E) G D with J > 0. 



Proof. Differentiating T with respect to E and J we obtain 

dT dyi dy2 dT dyi dy2 , . 



where 



/-/2 l + cos2(^) ^ ^ ^ 1 + Sin2(v9) , 

Jo (s(v?) - ys)^/^ Jo (•5(95) - 2/3)^/^ 



25 



E 




We remark that (/? i— > {s{(p) — 1/3) is strictly decreasing over [0,7r/2]. As 

Ai-A2 = V2 -j^dip, (4.7) 

Jo (5(99) - ya)-^/^ 



it follows from Lemma l2 . 51 that Ai > A2 > 0. On the other hand, since 
dyi 4yi dyi 4J 



1,2,3, (4.8) 



dE -3yf - 4y, + AE ' 8 J -3yf - 4yi + 4E ' 
it is not difficult to verify that 1)2. 12() still holds, whereas ()2.13|) is replaced by 

In particular, if J > 0, we obtain 

which proves ii). 

The main difficulty is the proof of i), which requires a more sophisticated argument. As yi — ys < 
s{ip) — ys < y2 — ys for all G (0, 7r/2), it is clear that 

< A A < ^ where c — '^'^^ 

{y2 - ys)^/^ ' (yi - ya)^/^ ' 4 

Using the upper bound for Ai and the lower bound for A2 we obtain from 1)4. 5|) : 

^ ^ + 2 ^ (49) 

(yi - y3)3/2 ai? ^ (y2 - ysf/^ BE' ^ ' ' 

Replacing the derivatives dyi/dE and dy2/dE with their expressions (|4.8() . we see that the right- 
hand side of (|4.9|) is positive if and only if 

^'i =^ yi(y2 - yif'\AE - ?,yl - 4^2) + y2(yi - yzf'\AE - Sy? _ 4yi) > . 
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Now, since AE — Sy^ — 4^2 < and AE — 3y^ — 4yi > 0, this inequality is equivalent to 

yi|4S - 3y| - 4y2| ^/yi-y3\3/2 n ^ ~ 2/3 ^ , 

— ttt; i : — r < i where < < 1 . 

y2{4:E - 3yf - 4yi) Vys-ys^ y2 - Vs 

Clearly, a stronger inequality is obtained if we replace the exponent 3/2 by 2 in the right-hand 
side. Thus it is sufficient to show that P2 > 0, where 

P2 = yi{y2 - y^fi^E - 3yl - Ay^) + y2{yi - y3?{^E - Syf - Ayi) > . 

To do that, we recall that yi, y2, ys are the roots of the cubic equation y^ + 2y'^ — AEy + 2J^ = 0. 
In particular, we have yi + y2 + ys = —2 and yiy2y3 = — 2J^. These relations allow to eliminate 
the variables E and ys from the expression of P2, which thus becomes a function of yi, y2 only. It 
is then convenient to set yi = y — z and y2 = y + z, where < z < y. After a straightforward 
algebra, we obtain the final expression 



P2 = 8z2((2 + 3y)2 + z2(3 + 4y)^ , 



which shows that P2 > 0- This concludes the proof. □ 

Using Proposition 14.11 it is now easy to verify that the Hamiltonian system associated to 1)4. 1[) 
is non-degenerate, i.e. the determinant A( J, E) defined in 1)2. 18(1 is always nonzero. 

Proposition 4.2 For all {J,E) £ D we have A{J,E) > 0. 

Proof. Differentiating the renormalized phase ^ with respect to J and E, we obtain as in 1)2. 19(1 

(4.10) 



5^ _ dyi dy2 _ r. j^yi , „ jdy2 

dE-^'-^dE^^'-^dE' dJ - ^'^JJ^^'^JJ-^'' 



where 



Bi ^ V2 - ^V^3 + Vsi^)-y3 ^ 

B2 Jo l(-y3)3/2(sM-y3)i/2 {V^3 + V'<v>)-y3y 



1 



^f—iZ^J.(,n\-„.\'i \ sin^ w / ^' 



(-y3)l/2(s((^) - y3)3/2 + ^ s{^) - y^f \ siu^i^) 

B3 = 2V2 / ^=^= ^=^^ = — - . 

-'0 V-ysV^iv) -y3{V^y3 + v^if) -ys) 

Thus replacing (|4.5|) . (|4.1()|) into (|2.18|) and proceeding as in the proof of Proposition ESI we obtain 

AiJ,E) = (A,B,-A,B2)iy2-yi)^^-B,^ . 

As dT/dE < by Proposition Ell it is thus sufficient to verify that Ai = A2B1 - A1B2 < 0. This 
inequality can be established using the same arguments as in the the defocusing case. Indeed, if 
we define (for i = 1, 2) Aj, Bj by TTTIh and Bj by (IT^ . we find Ai = 2(-y3)-'^A2 where 

A2 = A2B1 - A1B2 = A2B1 - A1B2 . 
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Now, we observe that i— > ct((/?) = (1 — 7/3 ^s((/?))^/^ is increasing over [0, 7r/2], because 1/3 < 0. 
Using Lemma 1231 we deduce that Ai > A2 > and B2 > Bi > 0, hence A2 < 0. □ 

As in Section |2 Propositions 14.11 14.21 ahow to determine the range of values of the period T 
and the renormahzed phase ^. We find that (T, : D ^ Z) is a smooth diffeomorphism, where 

Z> = |(T,^') e < T < vr, l^-l < ^(r)| , where ^{T) = -k - ^^'^ ^''^ . 

The domains D and D are represented in Fig. 5. 

Now, we fix (J, E) £ D and we study the stabihty of the periodic wave Uj^e{x, t) = e~'*VFj^£;(x), 
where Wj^e is a solution of (|4.1|) satisfying (|4.2j) . As in Section|31 we set Wj^e{x) = e'P^(5j,_E(2/cx), 
where k, i are defined in (|2.27|) and p = k + i. The discussion follows exactly the same lines as in 
the defocusing case, so we shall just mention the main differences. The function Q{z,t) defined in 
(|3.1|1 satisfies the evolution equation 

iQt + 4:ipkQ, + 4k'^Q,, + (1 - /)Q + |q|2q = , (4.11) 

and the corresponding energy functional reads 

£{Q) = £\2k'\Q,iz)\' - \\Qizt) dz . (4.12) 

In particular, if we define Sj^e by the second variation becomes 

Hj,E = £Ie^Qj,e) = -^k^d,, - Aipkd, - (l-p2) _ \Qj^e\^ - 2Qj^E ® Qj,e ■ 

As in the defocusing case, we rely on the result found for small waves in |12| Remark A. 2] and 
conclude that, when J = and > is sufficiently small, the operator Hj^e acting on Xq has 
exactly three eigenvalues {Ai, A2, 0} in a neighborhood of the origin, where 

Xi{E) = -3E + OiE"^) , \2{E) = E + 0{E^) , as E ^ . 

The other eigenvalues of Hj^e are positive and bounded away from zero. On the other hand, we 
know that zero is an eigenvalue of Hj^e of multiplicity exactly two for all {J,E) £ D. Indeed, 
this property was established in Section |31 by a general argument which uses only the symmetries 
of 1)2.1(1 and not the particular form of the nonlinearity. By continuity, it follows that Hj^e has 
exactly one negative eigenvalue for all {J,E) £ D, so that Proposition 13.21 remains valid in the 
present case. 

Finally, we recall that the equilibrium Qj^e of (|4.11|) is a member of a two-parameter family of 
travelling and rotating waves of the form 1(3. 11(1 . For sufficiently small {uj,c), the wave profile is 
given by Q^j'^{z) = XQji^e'{z), where A is defined by (|3.12|) and {J',E') £ D is the only point in 
a neighborhood of (J, E) such that (|3.13() holds. Let TCj^e be the Hessian matrix of the function 
dj,E defined by ()3.18|) . ()3.16|) . In view of Proposition 13. 2[ the general results of imply: 

Proposition 4.3 For all {J,E) £ D such that det(7ij^E) < 0, the periodic wave Qj^e is a stable 
equilibrium of J^. j j| ) in the sense of Proposition lA'. il 
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5 10 15 20 

Fig. 6: The quantity — det{Tij,E) is represented as a function of iJ G [E^{J), 20] for J = 0, 1, . . . , 5. Similar curves 
are obtained for higher values of J, thus indicating that det{Ti.j.E) is always negative. 

Proposition 14.31 is a conditional stability result, since it applies under the assumption that 
det{TCj^E) < 0. This condition is satisfied at least for {J,E) sufficiently close to zero, because we 
know from ^2] that 

nj,E = ^(J^ + , as {J,E)^ (0,0) mD . 

Moreover, the Hessian matrix TCj^e is easy to evaluate numerically for any {J,E) G D, and its 
determinant appears to be always negative (see Fig. 6). Thus we conjecture that det{Ti.j^E) < 
for all (J, E) G D. This property can probably be established rigorously using similar arguments 
as in the proof of Proposition 13. 5| but the modifications are not straightforward. 

4.2 Corotating waves (lo = —1) 

Finally we study the corotating waves of the focusing nonlinear Schrodinger equation. Our starting 
point is the stationary equation 

W^^{x) - W{x) + \W{x)\'^W{x) = , X G M . (4.13) 

The invariants of this Hamiltonian system have the following expressions: 

J = lm(WW,) , E = ^\W,\^ - ^\W\^ + ^\W\^ . (4.14) 

It is convenient to use the parametrization J = q{l + g^), where g G M. If J 7^ 0, the effective 
potential Vj(r) = J^/(2r^) — + r^/4 has a unique critical point at r = = y^l + q^, where 
Vj attains its global minimum: 

EM) = VjiVTT^) = i((?2 + l)(3q' - 1) . 

In that case, Ea. ()4.13|) has quasi-periodic solutions for all E > E^{J). If J = 0, the double- well 
potential Vo{r) has two minima at r = ±1 and a local maximum at r = 0. It follows that 1)4. 13|) 
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has (real) periodic solutions if— 1/4 < E < (dnoidal waves) and ii E > (cnoidal waves). 
Summarizing, the parameter domain where quasi-periodic solutions of ()4.13|1 exist is 

D = {{J,E) \ E > E^{J)}\ {{0,0)} . 

In the exceptional case {J,E) = (0,0), Eq. (BTTTni has the pulse-like solution W{x) = \/2/cosh(x) 
which corresponds to the solitary wave of the focusing NLS equation. 




Fig. 7: Existence domain for the corotating waves of the focusing NLS equation. The origin (J, E) = (0, 0) 
corresponds to the sohtary wave. The half-hne r_|_ (cnoidal waves) and the segment r_ (dnoidal waves) are the 
discontinuity lines of the phases $ and respectively. 

If (J, E) S D, the period T, the phase increment and the renormalized phase ^ of the quasi- 
periodic solutions are given by the formulas (|4.3|) . (|4.4|) . where s{ip) = yi cos^{ip) + y2 siv?{ip) and 
2/3 < < yi < 2/2 are the roots of the cubic polynomial P{y) = —y^ + + 4i?y — 2,P . It is 
important to realize that the phase $ is continuous on D, except on the half-line 

r+ = {(0,^;) gm2|^>o} = {{j,e) &D\yi = o} , 

(see Fig. 7). Similarly, the renormalized phase ^ is continuous on D except for a 27r-jump on the 
line segment 

r„ = {(0,^) gM^I -1/4<^<0} = {{J,E) GZ)|y3 = 0} . 

We shall thus use either <I> or ^' depending on the parameter region under consideration. 

Unlike in the previous cases, the period T is no longer a monotone function of the energy E. 
This is intuitively clear, as we expect that T +oo as {J,E) (0,0). In fact, one can prove: 

Lemma 4.4 The period T{J, E) satisfies 

dT dT 

— (0,^) > /or - 1/4 < ^ < , and —{0,E)<0forE>0. 
oE oE 



Proof. As in 1)4. 5() . we have 



dE ^ dE ^ dE 
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where the coefficients Ai, A2 are defined in ()4.6() . The only difference with the previous case is 
that the quantities y 1,112, 2/3 which appear in H4.6I) are the roots of a different polynomial. If J = 
and E' > 0, then = 0, 7/2 = 1 + \/l + 4:E, 1/3 = 1- Vl + 4£', hence 



If J = and -1/4 <E <d, then yi = 1 - Vl + 4^;, ^2 = 1 + Vl + ^E, 7/3 = 0, so that 

dT 2 

because Ai > A2 by □ 

On the other hand, the period T is still a monotone function of J if J > 0. As in ProDosition l4.11 
one can prove that 

5^' dT 

_(J,E) = _(,7,E) = --(J.E)>0, ifj>0. 

The important quantity is again 

/ dT d± \ / dT m \ 

A{J,E) = det[ g g (J,ii;) =det g g (J,ii;). (4.15) 
\ dj dj / \ a.j a.j / 

Here the first determinant in the right-hand side is meaningful if {J,E) ^ r_|_, and the second one 
if {J,E) ^ r_. 

Proposition 4.5 For all {J,E) e D we have A{J,E) > 0. 

Proof. Fix {J,E) G D. If {dT/dE){J,E) < 0, we know from Lemma Ol that {J,E) ^ r_. Thus, 
using the second determinant in H4.15|l and proceeding exactly as in Proposition 14.21 we obtain 
A{J,E) > 0. 

We now assume that {dT/dE){J,E) > 0, so that {J,E) ^ r+ by Lemma 031 Differentiating 
the expression of ^ in 1)4.31) with respect to E and J, we obtain 



dE = -^"^dE - ^"^dE ' dJ = -^'^^ - ^"^^ + ' 



where 



2cos^((^) l + cos^(y.) x 

= ^/o l(,.)VKv^)-.3 ^ ^(^)(-M-^3)3/0 ^^ = ^"+^^-' 

/--/y 2sin^((^) , l+sin^((^) 
^2 = V 2 / I — — ^=^^ + 1 d99 = B21 + -B22 , 

r- r^'^ 1 $ 

Jo s(v?) V ys 

Then, using (|4.5I) . (|4.16|) . and the first determinant in (|4.15|) . we find as in Proposition 12.61 
A{J,E) = {A,B2-A2B,){y2-y,)^^+B,^ . 
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As (dT/dE) > by assumption, it is sufficient to verify that Ai '^^^42^1 — A1B2 > 0. 
To do that, we observe that Ai = An + A12 where 

An = ^2^11 - ^1^21 , A12 = ^2^12 - ^li322 , 

and we prove separately that An > and A12 > 0. Both inequahties are easy consequences 
of Lemma Indeed, An > is equivalent to {Ai-A2)/{Ai+A2) < (;Bn-'S2i)/(i3n+i32i), or 
explicitly 

cos2((^) -sin2((^) cos2(99) - sm\ip) 

Jo {s{v)-y3r/' ^ ^ Jo s{^f{s{^)-y;Yl^ ^ _ 

7r/2 1^ r'^l'^ 'Y 



A stronger inequality (without the factor 3 in the denominator of the left-hand side) is obtained 
from Lemma [231 bv choosing I = [0,7r/2], f{if) = cos^ [ip) — siv? {(p) , g{ip) = {s{ip)—y^) / {s{(p))'^ and 

dfi = —z —-—^ TTTTT , where AI — ' 



(SM- 2/3)3/2 ' (,(^)_y3)3/2 • 

Indeed fj fgdu > {Jj f dfj.)(fj g dn) because f,g are strictly decreasing over I. Thus An > 0, and 
the same argument with g{^) = l/s{ip) shows that A12 > 0. □ 

To conclude this section, we fix (J, E) £ D and we study the stability of the periodic wave 
Uj^E{x,t) = e^^Wj^E{x), where Wj^e = e'P'^Qj^Ei'^kx) is a solution of (|I?T^ satisfying (|I?THl . 
Setting U{x,t) = e^^P^~^^^Q{2kx,t), we obtain from with 7 = 1 the evolution equation 

iQt + 4ipkQ, + ik'^Q,, - (1 + /)Q + |Q|2q = . (4.17) 
The wave profile Qj^e is a critical point of the modified energy 

SjMQ) = £iQ) + {l+P^)N{Q) - 4pkM{Q) , 
where £ is defined in H4.12() and N,M in H3.4I) . The second variation of £j^e at Qj^e is 
Hj,E = £Ie{Qj,e) = -^k^d,, - Aipkd, + (1 V) - \QjM^ - 2Qj,E ® Qj,e ■ 

Again, one can prove that Proposition 13.21 still holds in the present case. The fact that zero is 
always a double eigenvalue of Hj^e is established as in Section |31 and a direct calculation for small 
amplitude periodic waves (in a neighborhood of a plane wave) shows that Hje has exactly one 
negative eigenvalue. 

As in the previous cases, the equilibrium Qj^e of (|4.11|) is a member of a two-parameter family 
of travelling and rotating waves of the form 1)3. 11(1 . For sufficiently small {oj^c), the wave profile is 
given by Qy^{z) = XQj\E'iz), where A is defined by 1)3. 12(1 and (J', E') £ D is the only point in a 
neighborhood of (J, E) such that 

= (1 +/) - a2(1 +p'2) , c = AX'^k'p' - 4kp . 

If "^J.B denotes the Hessian matrix of the function dj^E defined by (|3.18|) . (|3.16|) . the results of fH^ 
imply: 
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Proposition 4.6 For all {J,E) £ D such that det{TCj^E) < 0, the periodic wave Qj^e is a stable 
equilibrium of J^. j7| j in the sense of Proposition 1^. il 

As in Section I^H we conjecture that det{TCj^E) < for all {J, E) G D. This inequality is true 
at least for small amplitude periodic waves (in a neighborhood of a plane wave), and numerical 
calculations indicate that it remains valid over the whole parameter domain D. In the particular 
case where (J, E) S r_ (dnoidal waves), the orbital stability with respect to periodic perturbations 
has been established in pp. Remark that Proposition 14.61 does apply in the case where {J,E) £ 
(cnoidal waves), which is not covered by the results of ^j, but as is explained in the introduction 
this is because we use in fact a more restricted class of perturbations. 
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